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ABSTRACT 


This document presents a baseline study for a generalized 
spacecraft attitude calibration/verification system. It 
can be used to define software specifications for three 
major functions required by a mission: the pre-launch 

parameter observability and data collection strategy 
study; the in-flight sensor calibration; and the post- 
calibration attitude accuracy verification. Analytical 
considerations are given for both single-axis and three- 
axis spacecrafts. The three-axis attitudes considered 
include the inertial-pointing attitudes, the reference- 
pointing attitudes, and attitudes undergoing specific 
maneuvers. The attitude sensors and hardware considered 
include the Earth horizon sensors, the plane-field Sun 
sensors, the coarse and fine two-axis digital Sun sensors, 
the three-a<is magnetometers, the fixed-head star track- 
ers, and the Inertial reference gyros. A review of the 
calibration/verification procedure currently performed 
on the Solar Maximum Mission (SMM) is also presented. 
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SECTION 1 - INTRODUCTION 


1.1 BACKGROUND AND MOTIVATIONS 

The spacecraft sensor calibration/verifcation support 
for attitude determination and control will become the 
primary ground support activity during the coming decade. 
The calibration support includes the determination of 
the sensor characteristics and alignments from which the 
sensors can be correctly modeled and the attitude can be 
accurately computed from the sensor measurements. The 
verification support consists of the comparison between 
the attitude computed on-board with that computed on the 
ground and the indication of requirements for new calib- 
ration updates. In the past, this calibration/verifica- 
tion function was performed by tailor made software sys- 
tems implemented to meet the mission requirements for 
specific missions. Due to the similarities in missions 
and sensors, many of these efforts were redundantly per- 
formed. Thus, the cost of such activity can be largely 
reduced if a baseline is available from which the mission 
software specifications can be defined for all mission 
requiring calibration/verification support. 

The purpose of this baseline study can be summarized as 
follows: 

* To generalize the in-flight calibration/ver- 
ification procedures to reduce redundant ef- 
forts for missions carrying similar sensors. 

* To help in developing general-purposed soft- 
wares to perform pre-launch attitude and bias 
observability studies so that 

(1) data acquisition schemes to optimize bias 
observabilities can be planned prior to 
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the mission and schedules into the mission 
time line. 

(2) Sensor configuration and orientations 
can be planned prior to the mission to op- 
timize the attitude and bias determination 
geometries . 

(3) Unnecessary efforts attempting to improve 
the attitude accuracy to below the geometric- 
al limit can be eliminated. 

(4) Unnecessary efforts attempting to deter- 
mine more than enough biases can be reduced. 

(5) Unnecessary efforts attempting to deter- 
mine biases from inadequate or improper data 
sets can be reduced. 

(6) Reliable estimates for attitude determin- 
ation accuracies can be obtained. 

* To simplify the calibration/verification 
software developments and operations by 
unifying the basic calibration algorithms 
and by sharing as many of the software mod- 
ules as possible among sensors and missions. 

1.2 DOCUMENT OVERVIEW 

To achieve the above goals, a review of calibration/ 
verification procedure currently being performed on the 
Solar Maximum Mission (SMM) is given in Section 2 to 
study its accuracy, possible enhancements, and feasibil- 
ity on other missions. Based on the SMM procedure, a 
generalized baseline is then developed from which future 
calibration and verification functions can be fashioned. 
This baseline is presented in Section 3 of the document. 
From this baseline, software specifications for the 
following three major functions required by each mission 


can be defined: (1) Pre-launch parameter observability 

and data collection strategy study; (2) In-flight sensor 
calibration; and (3) post-calibration verification. 
Finally, the analytical basis of the baseline and the 
equations required for each of the software modules in- 
cluded in the baseline are presented in Section 4. In 
these analytical considerations, both the inertial and 
the reference-pointing spacecrafts are included. The 
attitude sensors and hardwares considered include the 
Earth horizon sensors, the plane-field Sun sensors, the 
coarse and fine two-axis digital Sun, sensors, the fixed- 
head Star trackers, the inertial reference gyros, and 
the three-axis magnetometers. 


SECTION 2 - CURRENT SMM CALIBRATION/VERIFICATION 

PROCEDURE 

The SMM is the first of the Multimission Modular Space- 
craft (MMS) series. The scientific objective of the 
mission is the study of the complex solar flare phenom- 
ena. The spacecraft was launched into a hear-circular 
orbit on February 14 , 1980. The attitude determination 
and control hardware onboard the spacecraft includes two 
coarse Sun sensors (CSS), one fine Sun sensor (FSS) , two 
fine pointing Sun sensors (FPSS) , two fixed-head star 
trackers (FHST) , twt? three-axis magnetometers (TAM) , and 
one inertial reference unit (IRU) which contains three 
gyros, four reaction wheels, six magnetic torquing coils, 
and one onboard computer (GBC) . To meet the attitude 
determination and control accuracy requirements, sensor ' 
calibration/verification activities were performed by 
the SMM Attitude Ground Support System Software (AGSSS) . 
The AGSSS consists of six independent subsystems. Two 
of them are responsible for the calibration/verification 
functions: the SMM Attitude Determination Subsystem 

(SMM/ ADS) ^ ' 3 and the SMM Fine Pointing Sun Sensor Off- 
Null Calibration Subsystem (SMM/FOCS) 3 . The SMM/ ADS 
was used to calibrate the three-axis magnetometers, the 
fixed-head star trackers, and the inertial reference 
gyros; and the SMM/FOCS was used to perform the FPSS 
off-null calibrations. 

, 4 

These sensor calibration activities are reviewed in the 
following subsections. Four areas are addressed for 
each type of sensors The bias parameters defined, the 
calibration method used, the calibration results obtained, 
and a discussion of results. Some general comments on 
each of the calibration procedures are also presented to 
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serve as a guideline for future sensor calibration con- 
siderations. 

2.1 THREE-AXIS MAGNETOMETER 

2.1.1 Calibration Procedures and Results 

The TAM measures the geomagnetic field in the spacecraft 
body coordinate which is used in conjunction with one of 
the Sun sensor measurements to provide a coarse attitude 
determination by the Coarse Attitude Determination Sub- 
system (CADS) in SMM/ADS. The TAM measurements are also 
used for momentum unloading. Thus, reliable magnetometer 
calibration is essential for both coarse attitude deter- 
mination and spacecraft momentum management. 

The CADS is also responsible for the TAM calibration 
function. The calibration parameter defined in this 
function is a constant bias vector on the measured geo- 
magnetic field. The CA5JS uses an iterative least-squares 
technique to minimize the differences between the magni- 
tudes of measured geomagnetic field and that computed 
from a mathematical model. The following conclusions 
were obtained as a result of the SMM magnetometer calib- 
ration activities: 

(1) The magnetometer bias determined from a short inter- 
val can be totally unreliable. 

(2) The magnetometer bias can be accurately determined 
by processing either a full orbit of data or two 10 min- 
ute data passes separated by approximately 30 minutes. 

(3) Properly determined magnetometer biases remain rela- 
tively stable over long time periods. 

2.1.2 Discussions and Comments 

Some general comments and discussions on the SMM magnet- 
ometer calibration procedures and results are given in 
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the following; 

(1) In general, two types of systematic errors are like- 
ly to affect the TAM measurements: the sensor misalign- 

ments and uncompensated magnetic field contaminations. 

The former can cause an orientation error on the measured 
geomagnetic field while the latter causes a residual bias 
on the measurement. In the current SMM procedure# only 
the latter is modeled and determined. No attempt was 
made in calibrating the TAM sensor alignment parameters. 
This is reasonable so far as attitude determination at 
inertial attitude is concerned because the effect on 
data due to the sensor misalignment cannot be disting- 
uished from that due to attitude errors. In other words# 
the sensor misalignment parameters can be determined on- 
ly when the attitude is determined through other means. 

In this case# the accuracy of the sensor alignment is 
the same as that of the attitude# and the sensor mis- 
alignment parameters so determined can certainly not be 
used to further improve the attitude accuracy. However# 
a good TAM calibration can always improve the perfor- 
mance of the momentum unloading. Thus# it would be de- 
sirable to include the capability of determining the 
* misalignment parameters in TAM calibration procedure so 
far as the momentum management is concerned. In this 
case# the TAM alignment can be determined when the atti- 
tude is accurately determined by FPSS and FHST. The 
calibrated TAM can then be used to assist the momentum 
unloading process. 

In fact, the effects due to TAM misalignments and atti- 
tude errors are likely to be separated if proper slew 
maneuvers are performed. However, performing slew man- 
euvers during the early phase of the mission may not be 
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permissible operationwise. 

(2) The constant magnetometer bias vector can be separ- 
ated from the attitude error because the former changes 
the magnitude of the measured geomagnetic field while 
the latter doss not. Therefore, it is a good choice to 
use the magnitude of the geomagnetic field as the obser- 
vable in the least-square filter. This enables the de- 
termination of the magnetometer bias regardless of the 
accuracy of the attitude. However, as the spacecraft 
position changes along the orbit, the direction of the 
geomagnetic field B also changes. Depending on the in- 
stantaneous direction of the geomagnetic field, a con- 
stant bias can give very different effects on its magni- 
tude. For instance, a bias in 2 direction gives maximum 
effect on the field magnitude, B, if B is long 2 direc- 
tion; but gives essentially no effect on B if B is along 
St or £ direction. Thus, very different results in the 
bias vector b can be obtained if short data passes 
are used. On the other hand, when data from a long per- 
iod in orbit is used, the variation in B readily tells 
the direction and magnitude of the bias b. The more 
variation of B, the better for determining b. This ex- 
plains the first two conclusions resulting from the SMM 
magnetometer bias determination. In fact, B goes through 
approximately a cycle every half orbit period (about 48 
minutes). Therefore, a data pass with data uniformly 
distributed over a half-orbit period should be ideal for 
determining the magnetometer bias. 

Actually, the amount of variations ir the geomagnetic 
field directions within a half-orbit period differs 
strongly depending on the relative directions among the 
Earth pole, 2, the orbit pole, §, and the geomagnetic 
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dipole, ft. For a given set of orbit parameters, the 
directions of 2 and § are fixed in the inertial space 
while the vector ft rotates about 2 once a day with an 
angular separation of about 11.4 degrees from 2. Figure 

2.1 shows the approximate variation in B during a half- 
orbit period when a simple dipole model is assumed for 
the geomagnetic field. The two cases shown in Figure 

2 .1 correspond to the two extreme cases when ft is copla- 
ner with § and 2. As can be seen from the figure, the 
amount of variation in B reaches maximum when ft is in 
the direction of ft 2 which is coplaner with § and 2 but 
away from §. Thus, the data pass chosen f~~ the TAM 
bias determination will be optimal if the orbit time is 
chosen to satisfy this geometrical condition. 

2.2 FIXED-HEAD STAR TRACKER 

2.2.1 Calibration Procedures and Results 

The FHSTs measure the star directions in the spacecraft 
body coordinate. These measurements, together with the 
FPSS measurements, are used by the Fine Attitude Deter- 
mination Subsystem (FADS) of SMM/ADS to provide fine 
attitude determinations. Since the FPSS measurements 
can only determine the pitch and yaw angles when the 
spacecraft roll axis is pointing within 0.5 degrees about 
the Sun, the FHST measurements are responsible for the 
roll angle determination for all attitudes as well as 
the fine attitude determinations when isK spacecraft is 
slewed away from the Sun. Thus, reliable FHST calibra- 
tions are essential for accurate attitude determinations. 

The FHST calibration function is also performed by FADS, 
Two types of parameters were calibrated for SMM: The 

scale factors used to transform the measurements from 
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VARIATION OF 8 IN HALF-ORBIT PERIOD 



S - ORBIT POLE 
Z - EARTH POLE 

m - GEOMAGNETIC DIPOLE MOMENT 


Figure 2-1. Variations in Geomagnetic Field Directions 
During Half -Orbit Period. 
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raw counts to angles, and the sensor misalignment para- 
meters. The methods used to determine these parameters 
and the calibration results are described below. 

2. 2. 1.1 Calibration Methods for Scale Factor Determina- 
tion 

Two methods were used in determining the FHST scale fac- 
tors: the angular separation method and the empirical 

partial derivative method. In the angular separation 
method, the angular separations between stars observed 
by FHST were compared with that computed from the cata- 
log stars. The scale factors, which were assumed to be 
constants, were then determined by minimizing the differ 
ences between observed and computed angular separations. 
In the empirical partial derivative method, the errors 
in scale factors were determined, from the discrepancies 
between the attitudes determined from FHST and FPSS 
through partial derivative expressions. In these expres 
sions, the scale factor errors were assumed to be the 
only error source that contributed to the attitude dis- 
crepancies. 

2 . 2 . 1 . 2 Calibration Methods for Misalignment Determin - 
ation 

The FHST misalignment parameters were determined through 
an empirical partial derivative method similar to that 
used in the scale factor determination. In this method, 
the misalignment angles were computed from the discrep- 
ancies between the attitudes determined from FHST and 
FPSS through partial derivative expressions. Here, the. 
FHST misalignment parameters were assumed to be the on- 
ly error source that contributed to the attitude discrep 
ancies. 
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2. 2.1.3 Calibration Results 


The following results were obtained from the SMM FHST 
calibration activities: 

(1) The scale factor determination in general did not im- 
prove the FHST attitude determination accuracy regard- 
less of the method used. 

(2) The misalignment determination did improve the ag- 
reement between the attitudes determined from FHST and 
FPSS. 

(3) Calibration in roll direction was not possible from 
the current method because the reference roll angle 
could not be obtained from FPSS. 

(4) The uncertainty in determination (misalignment 
about ^ FHST f the FHST boresight) is much higher than the 
uncertainty in determination (misalignment about - 

2 FHST ) * 

2.2.2 Discussions and Comments 

In the following, some general comments and discussions 
on the SMM FHST calibration procedures are summarized. 

(1) Among the two methods used in the scale factor de- 
termination, the angular separation method should be a 
better method to use because the angular separations 
between stars are independent of both the attitude and 
the sensor misalignments. Thus, using the angular sep- 
arations between stars as observables, the scale factors 
can, in principle, be determined regardless of the atti- 
tude and sensor alignment accuracies. 

(2) The attitude accuracy was not improved by the scale 
factor determination suggests that the scale factor er- 
ror is probably insignificant comparing with other error 
sources. 
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(3) The FHST attitude and bias determination results ob- 
tained from single tracker may have 'high uncertainties due 
to geometrical limitations. The uncertainty is especially 
high for the orientation of the spacecraft about the FHST 
boresight, as demonstrated by the high uncertainty in the 
ot^ determination. This situation is analogous to that of 
FPSS where the orientation of the spacecraft about the FPSS 
borsight (the roll angle) is totally unde ter imined. However, 
if both star cameras are used, the attitude and bias obser-^ 
vabilities can be highly improved because of the increase 

in geometrical variations. 

(4) The current FHST calibration procedure depends on the 
attitude obtained from FPSS and hence provides no calibra- 
tion in the roll direction. In fact, it may be possible 

to calibrate FHST independent of FPSS if proper slew maneuvers 
are performed within the mission contraints. It is impor- 
tant to calibrate the roll alighment of FHST because this 
is the only sensor which provideis the roll information. 

2.3 INERTIAL REFERENCE GYROS 

2.3.1 Calibration Procedures and Results 

The gyro assembly in IRU consists of three two-axis gyros. 

The outputs from three of the six axes are used in angular 
velocity calculations to specify the three components of 
the spacecraft's angular velocity vector. This angular 
velocity vector is then used to propagate the spacecraft 
attitude. Reliable gyro calibration is important to cor- 
rectly measure the change in attitudes. 

The gyro calibration function is performed by. the Gyro 
Calibration Subsystem (CALIB) of SMM/AOS. A total of 
twelve calibration parameters are defined in CALIB. 

Nine of them form the 3x3 misalignment/scale factor cor- 
rection matrix and the remaining three parameters give 
the gyro drift rate bias vector. 
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The calibration method used in CALIB is a least-squares 
technique which minimizes the difference between the 
propagated attitudes derived from gyro measurements and 
the reference attitudes determined from FPSS and FHST. 

Slew maneuvers were performed to support the gyro calib- 
rations. During the SMM mission, the gyro calibration 
activity was repeated along with the FHST and FPSS sen- 
sor calibrations. The results showed that the gyro cal- 
ibration parameters so detemined depend heavily on the 
results of FPSS and FHST calibration. 

2.3.2 Discussions and Comments 

The strong dependence of the gyro calibration results 
on the FHST and FPSS calibrations indicates that the 
true gyro calibration parameters probably have not yet 
been reached. This is because the desired long maneuvers 
for the gyro calibration purpose were not achieved due 
to mission and operational constraints. Thus, the cur- 
rent gyro calibration results depend heavily on the at- 
titude determination accuracies at both ends of the cal- 
ibration slews. In fact, gyro measures the change in 
attitude which should be insensitive to the attitude un- 
certainties Caused by systematic errors as long as the 
same sensor was used in determining the starting and 
ending attitudes. Thus, the gyro may be better calibra- 
ted if the attitudes at both ends of a calibration slew 
are determined from the same sensor. However, it is impor- 
tant to collect as many observations as possible to reduce 
the effect due to random errors. 

2.4 FINE POINTING SUN SENSOR 

2.4.1 Calibration Procedures and Results 

The FPSS measures the projected Sun angles in the SMM 
xy and xz planes. These measurements are used in FADS 


2-10 


to determine the pitch and yaw angles of the attitude. 

To ensure that the accuracy of the FPSS readings over 
the instrument's field-of-view are consistent with the 
SMM mission requirements of ±5 arc-sec, the FPSS were 
calibrated in-flight by the FPSS Off-Null Calibration 
System CSMM/FOCS) . 

The FPSS calibration parameters consist of a set of 
eight transfer function coefficients for each of the two 
measured angles, at and A recursive least-square fil- 
ter was used in SMM/FOCS to determine these parameters 
by minimizing the residuals between the pitch and yaw 
angles computed from the FPSS measurements and the cor- 
responding reference angles obtained from gyro measure- 
ments for a set of slew maneuvers. 

The following conclusions resulted from the FPSS calib- 
ration activities. 

(1) The discrete slew data gave better results than the 
continuous slew data due to the preaveraging process 
perfromed on the latter which reduces the effect due to 
random errors. 

(2) The FPSS calibration results depend on the accuracy 
of the gyro calibration. 

* t 

(3) Good agreement between the attitude changes measured 
by gyro and by FPSS were obtained after the FPSS calib- 
ration. 

2.4.2 Discussions and Comments 

Since the FPSS calibration was based on the gyro measure- 
ments, the accuracy of the FPSS calibration depends 
heavily on the gyro accuracies. Good agreement between 
the attitude changes measured by gyro and by FPSS does 
not necessarily ensure good FPSS calibration. The cur- 
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rent FPSS calibration result* can be improved if the 
gyro can be better calibrated as discussed in Section 
2.3.2. Furthermore, many of the parameters attempted 
to be determined in the FPSS calibration are not really 
distinguishable. Actually, these parameters are obtained 
from power series expansions of some analytical expres- 
sions^ hence are not all independent parameters. 

They are actually functions of a few constants which 
carry certain physical meanings. Thus reducing the num- 
ber of parameters to be determined and defining proper 
relationships among the FPSS coefficients are essential 
in improving the FPSS calibration. This is given in 
Section 4. 
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SECTION 3 •» GENERALIZED CALIBRATION/VERIFICATION 

SYSTEM BASELINE 


This section presents the generalized system baseline 
for spacecraft calibration/verification functions. Sec- 
tion 3.1 describes the major functions which can be de- 
fined from this baseline. Section 3.2 gives an overview 
of the algorithm to be used in the baseline and presents 
the baseline in terms of functional baseline and block 
diagrams. The external interfaces between subsystems 
and sample data file structures are also provided. Sec- 
tion 3.3 then summarizes the software modules required 
by each of the functions and outlines the modules which 
are sharable by different functions and missions. 

3.1 MAJOR CALIBRATION/VERIFICATION FUNCTIONS 

To ensure accurate attitude determination and control, 
three major activities related to sensor calibration 
are generally required by each mission. Namely, the 
prelaunch attitude and bias observability and data col- 
lection strategy studies; the in-flight sensor calibra- 
tion and bias determinations; and the post-calibration 
attitude accuracy verifications. 

The first function requires a predictor which predicts 
the attitude and sensor bias determination accuracies 
achievable under various geometrical conditions. This 
function is important because the results of these 
studies often can be used to enhance the sensor selec- 
tions and configurations, to improve the mission time- 
line schedules, and to aid the data acquisition and 
maneuver plannings. Most of all, they provide the pre- 
launch knowledge of the attitude and sensor bias 
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determination accuracies achievable from each of the 
sensors at various geometrical conditions. This know- 
ledge is essential in ensuring the fullfillment of the 
attitude accuracy requirements, optimizing the attitude 
determination and sensor calibration techniques and 
procedures, simplifying the operational supports, and 
improving the understanding of the realtime attitude and 
sensor bias determination results. 

The second function requires a data processor which pro- 
cesses the realtime spacecraft data obtained from var- 
ious sensors to determine the spacecraft attitudes and 
sensor biases. This sensor calibration activity is 
usually necessary to achieve the mission requirements 
on attitude determination and control accuracies. With 
the help of the predictor, data can be acquired under 
the most favor. ->le conditions or from pre-scheduled 
maneuvers to optimize the bias observabilities. Since 
the sensor calibration function requires relatively long 
data spans, it is generally performed on ground with the 
results uplinked to the spacecraft. The onboard compu- 
ter then determines the spacecraft attitudes momentarily 
using the calibrated sensors. 

The third function requires a quality assurance system 
which compares the attitudes determined momentarily on- 
board with the attitudes obtained from the ground proces 
sing when sensors are recalibrated using long data spans 
This serves as a verification of the validities of the 
sensor calibration results carried onboard. If the dif- 
ference between the attitudes exceeds the attitude 
accuracy requirement, it indicates that the current sen- 
sor calibration information carried onboard needs up- 
dates and refinements. This verification function is 
important when the sensor biases are not truely constant 
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3.2 SYSTEM BASELINE 


This section presents the baseline of the generalized 
calibration/verification system (GCS) . From this base- 
line , the system specifications for the three functions 
described in Section 3.1 can be defined. An overview of 
the algorithm to be used by GCS is given in Section 3.2.1. 

The functional baseline and block diagrams are presented 
in Section 3.2.2. 

3.2.1 System Algorithm Overview 

The theoretical basis of GCS is a least-squares filter 
which determines the elements of a "state vector" by 
minimizing the square of the difference between the ob- 
served data and the expected data computed from an ob- 
servation model. Detailed descriptions, categorizations, 
and comparisons of least-squares filters can be found 
in Reference 5. 

The following least-squares estimators are chosen for 
the implementation of GCS. In the case of an inertially 
fixed state, a "batch least-squares estimator" will be 
used which updates the state vector after processing a 
block of observations and iterates the state until a 
converged solution is obtained. In the case of dynamic- 
ally varying states, a "Kalman filter" will be used to 
propagate the state from one time to the next. In a 
Kalman filter, the filter's confidence in its estimate 
of the state is allowed to degrade from one update to 

\ 

another using models of noise in the state vector. This j 

causes the influence of earlier data on the current 

state to fade with time so that the filter does not lose 1 

:3 

sensitivity to current observations. A Kalman filter 
reduces to a "recursive least-squares estimator" when 
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the state noise is zero. 

The mathematical formulations of the least-squares es- 
timator used in GCS is summarized in the following. 

To update the state and covariance during an inertial 
state, equations for a batch least-squares estimator 
are used. That is 


5 k 


5k-i + M 2 
^ 1 0 


(H T R - 1 v) 

c £-k-i 


-i 


p 0 <*o‘Sk-i>l 


(3-1) 


p = [p ^ 1 + z (h t r _ 1 h) ) _1 

K A-o ^-k-i 


(3-2) 


where 

k = iteration number 
X = state vector 
P = covariance matrix 

H * observation partial derivative matrix 
R = observation error matrix 
v = observation residual vector 
l = time index 


To propagate the state and covariance for a dynamically 
varying state, equations for a Kalman filter are used. 
That is, 


-i+i = *i-i ‘ (3-3) 

P l + . - *l P i^ + Wl (3-4) 

where 

i = time index 

L = state vector 

P = covariance matrix 

<P - state transition matrix 
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Q » state noise matrix 
P » state noise transition matrix 

The analytical expressions of the entries in Equations 
(3-1) to (3-4) for each of the observation and transi- 
tion models will be given in Section 4. 

Equations (3-1) to (3-4) are necessary for both calibra- 
tion and verification functions where data processing 
capability is required. However, for an observability 
prediction system, only Equations (3-2) and (3-4) are 
needed. In this case, a nominal state will be assumed 
which either stays constant or propagates with time 
through given functional variations. Furthermore, a 
predictor does not require observation data. It only 
requires an observation model which predicts the times 
and uncertainties of the observations. These differences 
due to different functions are indicated in the baseline 
whenever necessary. 

3.2.2 System Baseline and Block Diagrams 

Figure 3-1 shows the functional baseline diagram of 
GCS and its external interfaces. GCS will operate under 
the Graphics Executive Support System (GESS) and use the 
core allocation/deallocation graphics displays and 
interactive processing services provided by this execu- 
tive. The system contains four major parts: the main 

driver (DRIVER) , the state initialization and data prep- 
aration subsystem (SIDP) , the state and covariance up- 
date and propagation subsystem (FILTR) , and the attitude 
verification subsystem (VERIF) . The functional block 
diagrams of these four parts are shown in Figures (3-2) 
to (3-5). In these figures, software modules which can 
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Figure 3-1. Functional Baseline Diagram of GCS 
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be applied in general to all missions and sensors are en- 
closed by solid boxes and those modules which are mission 
or sensor dependent are enclosed with dashed boxes. 

The primary functions of the DRIVER is to handle inputs/ 
outputs and to serve as a linkage between the 6ESS exec- 
utive and the three subsystems of GCS. It reads the 
input parameters through either cards or graphic termin- 
als. The outputs can lead to the lineprinters, the 
graphic terminals , or the disc data files. Two disc data 
files can be updated by DRIVER as a result of the atti- 
tude ^nd bias determination: the Sensor Calibration 

Data Set which stores the updated sensor calibration 
parameter, and the Attitude Status File which stores the 
attitude determination results. 

The SIDP subsystem is responsible for the state and co- 
variance initializations and data preparations. To 
generalize the system, a State/Observation Definition and 
Linkage File (SODL) is required to define the observation 
and transition models included in the system and to link 
each model with its related state elements. A sample 
file structure for SODL is shown in Table 3-1. With SODL, 
the system will automatically define the corresponding 
state vector elements for the user once a given combin- 
ation of observation and transition models are selected. 

After the state vector is defined, and the initial co- 
variance and state are acquired, SIDP then creates a 
Standard Data File (SDF) from the individual raw data 
files or from given data modelings. The SDF has a stan- 
dard file structure which is required by the FILTR sub- 
system. Briefly, it combines various observation and 
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Figure 3-2. Functional Block Diagram of GCS/DRIVER 
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Figure 3-3. Functional Block Diagram of GCS/SIDP 
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Figure 3-5. Functional Block Diagram of GCS/VERIF 
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Table 3-1. Sample File Structure of SODL 


Record 

Number 

Record 

Description 

Word 

Number 

Word 

Description 

1 

Observation 
model header 

1 

Total number of observation models 

2-3 

Name of observation model 1 

4 

Record number of the information record for observation 1 

5 on 

Repeat words 2-4 for other observation models 

2 

Transition model 
header 

1 

Total number of transition models 

2-3 

Name of transition model 1 

4 

Record number of the information record for transition 1 

5 on 

Repeat words 2-4 for other transition models 

3 

' 

State element 
header 

1 

Total number of state elements 

2-3 

Name of state element 1 

4 on 

Repeat words 2-3 for other state elements 

4 

Information 
record for 
Observation 1 

1 

Total number of state elements related to this observation 

2 

State element ID of the first related element 

3 on 

Repeat word 2 for other related state elements 

5 

Repeat Record 4 
for other obs. 





















transition data into data blocks. Each data block is 
identified as either an observation block which measures 
an inertial state , or a transition block which measures 
the state transitions. A sample SOF structure is shown 
in Table 3-2. In the case of a predictor, the state ob- 
servabilities can be predicted either from actual data 
or from some data modelings. In the latter case, only 
the header information of SDF, which gives the time and 
the amount of data,, needs to be generated. There is no 
need of simulating the actual data content. 

The FILTR subsystem is the center of GCS,. It reads the 
data from SDF and performs all the necessary computations 
to optimize the state and determine the covariance. It 
processes one block of data at a time. Equations (3-1) 
and (3-2) are used to process the observation data blocks 
and Equations (3-3) and (3-4) are used to process the 
transition data blocks. The state is updated after each 
observation processing and propagated after each trans- 
ition processing. The batch least-squares filter of 
each observation block contains only a single time frame 
and only one iteration is performed on each data proces- 
sing. 

In the case of a predictor, the state will be kept at 
its nominal values and only the covariance matrix needs 
to be updated and propagated. 

The VERIF subsystem is required by the verification func- 
tion only. It reads in the attitudes to be verified from 
the attitude status file, compares these attitudes with 
the attitudes obtained from FILTR when the updated sensor 
calibration parameters are used and sends an indication 
flag when the difference is greater than the attitude 
accuracy requirement. 
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3.3 SUMMARY OF SOFTWARE MODULES 


A summary of functional software modules required by GCS 
is given in Table 3-3. For each of the modules, Table 
3-3 specifies the functions by which the module is 
required and indicates if the module is mission or sen- 
sor dependent. This provides a guideline of the amount 
of software sharable by different missions and functions 
if the generalized calibration/verification system is 
implemented. 


Table 3-2. 


1 


u> 

l 


Record 

Number 


1 


2 


3 


Record 

Description 


File Header 


Header of 1st data block 


Header of 1st time frame 
in the block 


4 to n Data records of 1st time 
frame in the block 


n+1 to n Repeat records 3 to n for 
other time frames in the 
block 


nrt-1 
to end 


Repeat records 2 to m for 
other data blocks 


ample Fxle Structure of SDF 


Word 

Number 


Word 

Description 


Total number of data blocks 


File Type (1= headers + data, 
2= headers only) 


Block ID (1= observation, 
2= transition) 


Number of time frames in the block 


1-2 


Frame time 


Number of data record in the frame 


5 on 


Number of data points in model 1 
Repeat word 4 for other models 











Table 3-3. Summary of Functional Software Modules 


Required by 


Functional 

Subsystem Software Module 


SESS and other 

external 

interfaces 


I/O of general 


Pred- 

iction 



State element 
determination 
from SODL 


State and 

covariance 

initialization 


SDF 

generalization 
from raw data 
files 


SDF header 
generation from 
data modeling 



VERIF 


tate and 
covariance update 
and propagation 


Observation error 
and partial 
derivative 
matrices compu- 
tations 


Observation 
residual compu- 
tation 


btate and covar- 
iance transition 
matrices compu- 
tations 


ASF read 


Attitude 

computations 




X 

No 

X 

No 





































































SECTION 4 - ANALYTICAL CONSIDERATIONS 


This section presents the analytical equations required 
by the Generalized Calibration/Verification System (GCS) 
described in Section 3. Three types of attitudes are 
considered: (1) the single-axis spin-stabilized space- 

craft , (2) the inertially-fixed three-axis spacecraft, 
and (3) the reference-pointing three-axis spacecraft. 

In the case of the single-axis spin-stabilized space- 
craft, the attitude is defined by the orientation of the 
spin-axis which is inertially fixed. To determine such 
an attitude, the plane-field Sun sensor which measures 
the Sun angle, and the Earth horizon sensor which meas- 
ures the time between the Sun-sighting event and the 
Earth horizon-crossing event are usually used. In the 
case of the three-axis stabilized spacecraft, the atti- 
tude is defined by the orientation between the space- 
craft body axes and some reference axes. The sensors 
required for three-axis attitude determination typically 
include the coarse and fine two-axis digital Sun sensors, 
the three-axis magnetometers, the fixed-head Star track- 
ers, and the inertial reference gyros. 

Section 4.1 gives a general description of the coordin- 
ate system and attitude definitions to be used in the 
section. The transformation between different attitude 
representations, the modeling of attitude maneuvers, and 
the computation of spacecraft spin rate for reference- 
pointing spacecrafts are also given in Section 4.1. Sec- 
tion 4.2 provides the analytical considerations for each 
of the seven sensor types mentioned above. For each 
sensor type, the following subjects will be covered to 
include the necessary equations required by GCS: the 
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observation modeling, the state element definition, the 
nominal state, and the observation partial derivatives 
with respect to the state elements. For the case of gy- 
ros, equations required to perform state transitions and 
covariance propagations are also included. Section 4.3 
discusses the estimation of observation errors and the 
data modeling for a prediction system. Section 4.4.. 
describes the criterium used in attitude verification. 

4.1 COORDINATE SYSTEMS AND ATTITUDE DEFINITIONS 

4.1.1 Coordinate Systems 

For the purpos ? of spacecraft attitude determination and 
sensor calibration, four types of coordinate systems are 
frequently used. These are: (1) the geocentric iner- 

tial (GCI) coordinate, (2) the body-fixed coordinate, 

(3) the reference coordinate, and (4) the sensor coordin- 
ate. 

The GCI coordinate is fixed in the inertial space with 
x-axis pointing to the vernal equinox and the z-axis 
pointing to the celestial north pole. The body-fixed 
coordinate is a coordinate system which is fixed in the 
spacecraft body reference frame. The reference coordin- 
ate is needed for reference-pointing three-axis space- 
crafts only. It is a coordinate system defined by the 
orientation of the reference body. For example, the 
axes of the reference coordinate, of an Earth- 

pointing spacecraft expressed in the GCI coordinate is 
defined by 

0 = R/ j R | = local vertical 

0 = # x V/|0 x V| (4 _ 1} 

^ X ft/|0 X 

wh« 5 re R and V are the spacecraft position and velocity 
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vectors respectively. All vectors i,t\, Eq. (4-1) are ex- 
pressed in the GC1 coordinate. The reference coordinate 
for a Sun-pointing spacecraft such as SMM is defined by 
the SUN-coordinate as given in Reference 1. In general, 
the reference coordinate is not fixed inertially. The 
time variation of the reference frame is obtained from 
the ephemeris information. The sensor coordinate is 
defined by the sensor hardware mounted on the spacecraft. 
The orientation between the sensor coordinate and the 
body coordinate is usually described by an alignment 
matrix. Refinement of this alignment matrix is one of 
the major objectives of the in-flight spacecraft sensor 
calibrations. 

In the rest of the document, a vector expressed in these 
four coordinate systems will be subscribed by I, B, R, 
and S respectively. Vectors without Subscript are as- 
sumed to be expressed in the GCI coordinate. 

4.1.2 Attitude Definitions 

The definitions of attitudes for the single-axis and 
three-axis spacecrafts are given in the following sub- 
sections. 

4. 1.2.1 Single-Axis Attitude 

The single-axis attitude is defined by the right ascen- 
sion and declination of the spacecraft spin axis in the 
GCI coordinate. That is, 

A * cos a cos { Xj + sin a cos 5 Yj + sin <S (4-2) 

where 

A * unit vector along the spin-axis 

X-, Y_, Zj are the three unit vectors in GCI 
coordinate 

a * the right ascension of attitude 
6 * the declination of attitude 
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4. 1.2.2 Three-Axis Attitude 


The three-axis attitude is in general defined by the or- 
ientation of the spacecraft body coordinate relative to 
the reference coordinate. This orientation can be repre- 
sented by a coordinate transformation from the reference 
coordinate to the body coordinate as given below. 

— B * I B 1 £r <4-3) 

where r fi and r R corresponds to any vector y represented 
in body and reference coordinates respectively. The 
matrix [b] which defines the attitude is generally ex- 

. 5 

pressed m terms of Euler angles. Different conventions 
have been used in defining the Euler angle representa- 
tions of attitudes. To standardize the attitude defini- 
tion for GCS, the roll (<f>) , pitch (0), and yaw (i|>) angles 
of an attitude are defined as the 1-2-3 Euler rotation 
angles of B . That is, 



The Euler angles can be determined from the components 
of [b] by the following equation. 

<J> = tan -1 ( “® 3 2 / B 33 ) 0<<)><2tt 

0 = sin” 1 (B 31 ) -£<0<f (4-5) 

= tan” 1 (“ b 2 i/ b h) 0<\p<2it 

Another frequently used representation for three-axis 
attitude is the quaternion representation. 5 The quater- 
nion representation does not involve trigonometric func- 
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tions and therefore is convenient to use in computations. 
Also because it carries simple product rule for succes- 
sive rotations, the quaternion representation is usually 
used in computing the attitude transitions. In most 
practice, the quaternion representation of attitude is 
used to express the orientation of the spacecraft body 
coordinate relative to the GCI coordinate. That is, 


= [BjlA]^ (4-6) 

= [c(q)]r x 


where (b) = the transformation matrix from reference 
coordinate to body coordinate as given in 


Eq. (4-4) 

[A] = the transformation matrix from the GCI 
coordinate to the reference coordinate 

[Cj = the transformation matrix from the GCI 
coordinate to the body coordinate 

_ m t 

<3 = [q,q 4 ] is the quaternion corresponding to 

the matrix ( C ] . 

The matrix [A] can be obtained from the ephemeris infor- 
mation by the following: 


[A] = (I) 


for inertial-pointing spacecraft 



(4-7) 

for reference-pointing spacecraft 
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where Uj, v Wj are the three unit axis vectors of the 
reference coordinate expressed in the 6CI coordinate. 


The expression of f C] 

in terms of the quaternion is 

en by 

r q*-q*-q*+q* 

1 2 3 4 

2 (q q +q q ) 

2 (q q -q q ) 


12 3 4 

1 3 2 4 

fc (q)l = 

2 (q q -q q ) 

- q 2+q 2 -q2+q 2 

2(q q +q,q ) 

1 1 

12 3 4 

12 3 4 

2 3 14 


2(q,q +q,q L ) 

13 2 4 

2 (q q ~q q ) 

2 3 14 

-q 2 -q 2 +q 2 +q 2 

1 2 3 4 J 


To determine the quaternion components in terms of the 
elements of (c), we have 


5, “ ‘»l 1+C H +C 22 +C 33l 

q = 1 (C -C ] 

3 4q 4 lz Zl 

q = 1 [C -C 1 

2 4q 4 31 13 

q = 1 [ C -C ] 

l 23 32 


where 


(4-9a) 


q 2 +q 2 +q 2 +q 2 = 1 

1 2 3 4 


(4-9b) 


The quaternions carry a convenient product rule for suc- 
cessive rotations. Namely, if 


[C(q " )] = [ C (q )] [C(i)j 


(4-10a) 


then 
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q 


ii 




(4-10b) 


If both EuJ er angle and quaternion representations are used 
for the same transformation matrix, the Euler angles can 
be obtained directly from the quaternion components by 


4> = tan ^ 

'-2q 2 q 3 +2qi.q4 

, -q 2 -q 2 +q 2 +q 2 


1 1 2 3 4 

n . -1 
0 = SIR 

(2q q +2q q ) 

13 2 4 

, ^ -1 
^ = tan 

f-2qiq 2 +2q3q4 

1 

l 


0 < ij> < 2ir 

-I s - 6 s 1 

0 £ ¥ < 2it 
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4.2 ANALYTICAL CONSIDERATIONS FOR ATTITUDE SENSORS 

The analytical equations required by GCS are given in the 
following subsections for each of the seven sensor types 
of concern, namely, the Earth horizon sensor, the plane- 
field Sun sensor, the coarse two-axis digital Sun sensor, 
the fine two-avis digital Sun sensor, the three-axis mag- 
netometer, the fixed-head star tracker, and the inertial 
reference gyro. The first two sensor types are applied 
to the single-axis spacecrafts, while the remaining five 
sensor types are used for the three-axis spacecrafts. 
Detailed descriptions of these sensors and their model- 
ings can be found in Reference 5. 
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4,2.1 Earth horizon Scanners 


The Earth horizon scanners measure the time durations 
between the Sun-sighting event and Earth horizon cross- 
ing events for the spinning spacecrafts. The equations 
for the sensor modeling and partial derivatives are sum- 
marized in the following. Detailed derivation of equa- 
tions for the Earth horizon scanners can be found in 
Reference 6. 

The analytical expression for the horizon-in/out time 
model is given by 

v = — (<t> - <|>„ “ A<|>- A6 + A <p a +2irn) n=0,+l (4-12a) 

■* i,\ rl S n 


where 

0) 

<t> 


spacecraft spin rate 

Sun to horizon-in/out dihedral angle 


= tan 


-1 


(S x H) 


A A A A A 

H - (S • A) (H * A) 


(4— 12b) 


A = spin axis attitude as given by Eq. (4-2) 

A 

S = Sun unit vector 

A 

H = unit vector along horizon line of sight at hori 


zon- in/out time 

A a a 

.= cos (p+Ap) E + sin (p+Ap) (M sin A+ N cosA) (4-12c) 

A 

E = nadir vector 
M = A x E/sin n (n=nadir angle) 

/\ A A 

N = E x M 

cos A = cos (V+Ay ) -cos (p+Ap) cos n { . + . for horizon-oat, 
sin (p+Ap) sin n for horizon-in) 


where y=sensor mounting angle 

Ay=sensor moutning angle bias 
p=angular radius of Earth 
Ap= bias on angular radius of Earth 


if 
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sin A * ± (1-cos 2 A) 55 (' + ' for horizon-out, ' for in) 

<j> * azimuth mounting angle between Sun sensor and 

horizon sensor 

Aij) - constant azimuth bias due to mounting error or 
electronic delay 

A<J> S * azimuth error due to Sun sensor plane tilt 

= sin” 1 (cot 6 tan e_) (4-12d) 

5 

where 8=Sun angle 

s^sun sensor plane tilt 

= sin” 1 (cot (Y+Ay) tan e H ) (4-12e) 

where e =Earth sensor plane tilt 

n 

Thus, the state vector for the Earth horizon-in/out model 
can be defined as the following. 

X = (Yf <s, A4>, Ay# Ap, Eg, e H , to , At ) T i ~'' 

where At is the orbit in-track time error. 

The equations for the partial derivatives of the obser- 
vation with respect to the state elements are summarized 
in the following: 


lx 

3a 

lx 

35 


1 

' 3a 


to 


I 


30 


- secA(J> esc 3 tane S. ) 

0 s 3a 

A 

3 A ) 

- sec A<p a esc 0 tane s s * 


3y _ _1^ 

3A(j) ” to 

Iay = i III?- SeC A *H sine H Sin(Y+AY) ] 

lx _ 1 34 > 

3 p 0) 3 Ap 

1^- = -1 secA<t> 0 cot8 sec 2 e c 
3e„ oo S o 


4-9 


WrattaKwr*" 


lx m 

8e H 

1 

0) 

(it 

( S£ h 



- tan 

lx = 

-X 


3 co 

W 


lx = 

1 

d± 

3At 

0) 

3At 


'H 


'H 


i 


where 


84> _ 
3x 


= -L- .... / 1 12J. _9jl 9Qz\ 

1 l+(§‘) \ q 2 3x i 1 2 3 *i) 


= A • (S X H) 


A A 


Q 2 = S • H - (S • A) (H A) 


A A A A r\.j 

3 Q i = (S X H) • + (A x S) • 

3 *i i ^ *■ 

3Q a _ £ 3H " ■ 3 a . 

3x7 ^i. ’ ’ ' 

A A • 

+ (s • A) Jsin ( Y+Ay)cos e „ -3Aat 
l H jx. 

+ cos ( Y+ Ay) sin e H 

3A = X - I 1 Z (if X. = <*) 

3x I Aj I X 

-*!-%*! (if *i =6) 


= 0 (otherwise) 

1 (if Xj-AY) 

= 0 (otherwise) 


3 e H 


^ 1 U* «*- «„> 

= 0 (otherwise) 
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The partial derivatives of the horizon crossing vector 
with respect to the state elements, 3H/9x^, can be ob- 
tained by the following. 



9 H = H x E sin (Y+Ay ) cos e 
9A Y E • (H x A) 


iE = H x A sin (P+Ap ) 
9A P E • (A x H) 


9 H = H X E cos (Y+Ay ) sin e 

3 - A A A ** 

e H E • (H x A) 

4 ■ Ltl f(H - 6 IgfAp, . jl] 

9At E • (A x H) *• J 

where V is the spacecraft velocity. 

Two other observation models are oftenly used for the 
Earth horizon scanners, namely, the "Earth-width" model 
(y ) and the "Sun to Earth mid-scan time" model (y ) . 

The former models the time span between the horizon-in 
and horizon-out events and the latter models the average 
of the Sun to horizon-in time and the Sun to horizon-out 
time. The partial derivatives for these two models are 
given by 
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+ 


®*M - H( 
**i 




) 


where y ^ and y Q are the Sun to horizon-in time and Sun to 
horizon-out time models respectively and x i are the state 
elements . 

4.2.2 Plane-Field Sun Sensor 

The plane-field Sun sensor measures the Sun angle ($) be- 
tween the Sun-vector and the spacecraft spin axis. It 
is usually used in conjunction with the Earth horizon 
sensors to determine the spacecraft attitudes for single- 
axis spacecrafts. The sensor modeling and partial deriv- 
atives are summarized below. 

The Sun angle model for the plane-field Sun sensor is 
defined by the following 

y - cos -1 ( cos 3 sec e_) - A3 (4-14) 

where 3 is the Sun angle, e s is the Sun sensor plane 
tilt, and A 3 the gun angle bias. The state vector 
for this observation model is thus defined by 

x = (a, 6, A3, e s ) T (4-15) 

The partial derivatives of the observation with respect 
to the state elements are 

ay = x(a 2 si- ais 2 ) 

3a 

Sy = X (cos3 As- S 3 ) 
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where 

X - 1 

^coa z e g -coa 2 0 


4.2.3 Coarse Two - Ax i 3 Digital Sun Sensor 

The coarse two**axis digital Sun sensor measures the pro- 
jected Sun angles, oc and 3, in the xy and xz planes of 
the sensor coordinate. The measurements are related to 
these angles by the following: 


y ct * N a * |< a+b a ) (4-16a) 

* N 6 * lg (0+b 6 ) (4-16b) 

where K a and are the scale factors and b a and bg are 
the Sun agnle biases. Thea and 3 angles are related to 
the Sun vector by 


a 


tan 


tan 


•1 s ya 
XS 

(4-17a) 

Vzs) 

(4-17b) 

S xs 


and S zs are tbe components of 

the unit 


Sun vector in the Sun sensor coordinate. To express S, 
in terms of the attitude, we have 



(4-18) 


s s - Im](b)s r 

where [ M ] is the Sun sensor alignment matrix which trans- 
forms a vector from the spacecraft body coordinate to the 
sensor coordinate , [b] is a function of the Euler angles 
which defines the spacecraft attitude (Eq.(4-4)), and § R 
is the Sun vector expressed in the reference coordinate 
which can be obtained from the ephemeris information us- 
ing Equations (4-6) and (4-7). Thus, the Sun sensor 
measurements N a and Ng are related to the Euler angles 
$, 6 , *P through the Equations (4-16) to (4-18) and (4-4). 

Assuming errors due to sensor misalignment are negligible, 
then the state vector for coarse two-axis Sun sensors can 
be defined as 


X — (4> , 9, ip , K^, b^) 


i=a, S 


(4-19) 


The partial derivatives of the observation with respect 
to the state Elements are given by 

i . Yi <i=ot, 3) 

iiq - jq 


3y. _ 1 
J i — n — 

3b. K i 

3y . _ 1 ,3a v 

^ t : ( ^ t ’ 


(i =a, 0) 


(*j = <j>, 0, i|»,) 


& =4? ( -i— > <*j ■ e - ♦'» 


3x . 
D 


where 


3a 




S 2 +S 2 
xs ys 




(4-20a) 
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1 


(*l-20b) 


36 

3x. 


xs+ zs 


h, 


and the partial derivatives of S with respect to the 
Euler angles are 


3S 3 
3x 3 

t«i 3 ;■> 
3x3 

- 

/V 

S R 





(4- 

■21) 

where 
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0 

-sin^ 
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0 

-sin<}> 

cos4> 
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sin0 

0 

COS0 
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-sin4> 
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: sin6 
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0 
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0 
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0 

-sin4> 
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L sin^ 

cosip 
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0 

-sin(? 


\ 

0 

0 s 

3 [B] .. 

"COSlp- 

■sini^ 
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0 

1 
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0 

COS<J> 

sintf) 

3 i|> 

0 

0 
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sin© 

0 

COS0 


0 

-sin<}> 

COS(j) 


4.2.4 Fine Two-axis Digital Sun Sensor 


Similar to the coarse Sun sensor, the fine Sun sensor al- 
so measures the a and 6 angles of the Sun vector. How- 
ever, more parameters are included in the observation 
models to describe fine residuals in the measurement. 

As shown in Reference 5, the fine Sun sensor measurements 
can be modeled by the following expressions 


L a 




tan a 
tan 6 


A1-A3 sin(k a tana) 
B1-B3 sin(kg tang) 


( 4-23a) 
(4-23b) 


where A,, B,, k a , and k s are the eeneor parameters to be 
calibrated. It can be shown that the following expres- 
sions for « and Sin terms of N a and N g are obtained when 
proper iterative approximations are performed. That is, 


a - tan” 1 A l +A 2 N a +X3Sin(AHN a +A5)+A«sin(A 7 N a +A.) 

i 

0 . tan" 1 B 1 +B 2 N g +B S 8in(B.,N 6 +Bs)+B*ain(B7N g +B.) 


(4-24a) 

(4-24b) 


In Equation (4-23), the high order coefficients are not 
independent parameters. They are related to the param- 
eters in Eq. (4-22) by 
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k a A 2 

B 4 “ 

k g B 2 
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The state vectors for the fine Sun sensor are defined as 


^3 * (<{> f 9 ,!|>,Ai ,A 2 ,A3 ,k a ) 

*0 a (<J) ,0 ,lj>,Bl ,B 2 ,B3 ,k g ) X 

The partial derivatives of y a are given below 
expressions can be obtained for that of y g . 



3Ai A 2 


«a 

3A 2 A 2 

3Y a — ^ — sin (^a tan a ) 

3A 3 2 


(4-26a) 

(4-26b) 

Analogous 
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a 
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9 Y a 
a X^ 


A 

cos (k a tana) tan a 

2 

1 

1 - k a A 3 cos (katana) aa 

A 2 cos z ct ^ 

9x3 


where Xj » <p,0,ip and 9a/9x3 is given by Equations (4-20) 
to (4-22) . 


4.2.5 Three-Axis Magnetometer 

The three-axis magnetometer measures the components of 

the geomagnetic field in the magnetometer frame, H_. 

— — s 

The sensor observation can be modeled as 


Z = H s + b 

= CM] [B] [A] Hj + b (4-27) 

where [M] is the magnetometer alignment matrix which trans- 
forms a vector from body coordinate to the magnetometer 
sensor coordinate, [B] and [A] are defined in Equations 
(4-4) and (4-7) respectively, Hj is the geomagnetic field 
in the GCI coordinate which can be computed from mathema- 
tical models of the geomagnetic field, and b is a constant 
bias vector on the magnetometer measurements. 

Assuming the error due to sensor misalignment is negligible, 
then the state vector for the magnetometer observations can 
be defined as 

X = ($,e,ifi, b) T (4-28) 

The observation partial derivatives with respect to the 
state elements are given by 
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(i»l,2,3) 


9 *i - 

■5B3” 



<5 


ij 


(M) 1UL [AJ H x 
3x, 


where 3[B|/3*i is giv**’ 'n Eq. (4-22). 


4.2.6 Fixed-Head Star Tracker 


Analogous to the two-axis Sun sensors , the star tracker 
measures the direction of a given star in the tracker 
coordinate. The two measurements are related to the 
star direction by the following equations 


yH 




H 

_ 1 _ 

3 V 


sm 


tan 


Zs 

F 


XL 

i* 

xs 


(4-29a) 

(4-29b) 
T 


where S„ and S„ are the scale factors and F * (F ,F ,F ) 
H v s xs ys zs 

is the unit star vector in the tracker coordinate. 


F g can then be expressed in terms of the attitude para- 
meters by 


F S " [M][B][A] * z (4-30) 


Here, (mJ is the star tracker alignment matrix which 
transforms a vector from the body coordinate to the 
tracker coordinate, (b] and f A) are defined in Equations 

A 

(4-4) and (4-7) , and F^ i& che unit star vector in the 
GCI coordinate which can be obtained from the star cata- 
log. In order to calibrate the star tracker alignment, 
the matrix M can be expressed in terms of three rota- 
tion angles as following. 


[M] = T a (6) T 2 (8) Tx (a) 

where , Tj, T 3 are the rotation matrices defined in 
Equation (4-4) . Thus, the state vector for the star 
tracker can be defined as 
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X = (<{>,e,!|>,a,B,$,s H ,s v ) (4-31) 

The partial derivative equations for the star trackers are: 
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where = ip , 6 , \p ,a, B f <S . The partial derivatives of F g with 
respect to x^ are given by 

= B H A 1 Fj (if X £ « <x,&,6) 


3F 


Txf = ^ M 1 [&] P x 


(if x^ = <p,e,ip) 


where 3 [B]/3x^ is given by Eq. (4-22) . Similar Equations 
can be derived for 3[M]/3x i . That is, 


urn - 

3 a 

cos 6 
-sin6 

sin60 

cos60 

cos 6 
0 

0 

1 

-sin8 

0 

. 0 

0 

1 

sin3 

0 

cos 8 

3 [M] = 
TT~ 

cos 6 

sin<50 

-sin8 

0 

-cos 8 

-sin6 

cos 60 

0 

0 

0 

0 

0 

1 

cos 8 

0 

-sin8 

3 [M] 

«> £ 

-sin6 

cos60 

cos 8 

0 

-sin8 

-cos 6 

sinoO 

0 

1 

0 

o 0 

0 

0 

0 

sin8 

0 

cos 8 


0 0 0 
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0 -cosa-sina 

10 0 
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0 -sina cosa 

10 0 ' 
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4.2.7 Inertial Reference Gyro 

The inertial reference gyros measure the spacecraft angular 
velocity relative to the inertial coordinate. This measure- 
ment can be used to compute the state transition and to pro- 
pagate the covariance matrix. It can also be used as an 
observation to calibrate the gyros. These two functions are 
separately in the following subsections. 

4.2. 7.1 Gvro Calibration 

The gyro measurements are related to the spacecraft an- 
gular velocity by the following relation: 

Y = [g] <w + b) (4-32) 

where [g] is the inverse of the gyro scale factor cor- 
rection and alignment foatrix, w is the averaged space- 
craft angular velocity over the gyro measurement frame, 
and b is the gyro drift rate bias. The angular velocity 
w can be determined from the initial and final attitudes 
using the quaternion representation. If q(t-t) and q(t) 
are the initial and final attitudes, then w is obtained 
from the following procedures. 

Let Q = q (t-t) * q(t) * (Q, (4-33a) 


where q* * the conjugate of q 

* (-q» q 4 ) T 


then, 

w 

where 

0 ) 


(D Q 

T 

sin (-j ) 

— cos” 1 Q. 
t 4 


(4-33b) 


(4-33c) 


Equation (4-33a) requires the quaternion multiplication 
which is given by Eq. (4-10b) . For the case of observ- 
ability predictions, the nominal spacecraft angular vel 
ocity can be obtained directly from the ephemeris in- 
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formation. The angular velocity for an Earth-pointing space- 
craft is given by 


V . u 



where R is the distance between the spacecraft and the 

* A 

Earth, V is the spacecraft velocity, and u, v are defined 
in Equation (4-1) . 


To model the attitude maneuvers starting from time to, 
the attitude at time t is given by 


[c<5 t ,] = [» t ] [aJ 

where A^. is given by the ephemeris and 
Equation (4-4) with 


(4-35) 

is given by 


4> = $ 0 +<H 

0 = 0+0 1 

® • 

tp = t 


for roll-maneuvers 
for pitch-maneuvers 
for yaw-maneuvers 


(4-36) 


where *^0 Q are the spacecraft roll-pitch-yaw angles 
at time t Q . With the quaternions determined from Equa- 
tions (4-35) and (4-36) , the angular velocity w corres- 
ponding to a given maneuver can then be obtained from 
Equation (4-33) . 

Assuming the errors in the initial and final attitudes 
give neglibible effect on the angular velocity computa- 
tion, then the state vector for the gyro measurement can 
be defined as 


X - ( G ij> (i = 1 to 3, j = 1 to 3 ) (4-37) 

where G i; . is the component of the matrix [G]and b. is the 
component of the vector b. 
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The partial derivative equations for the gyro measure- 
ments are: 

9 y i - (w. + b. ) (ifl to 3, j-1 to 3, 

' - 1 K * k*l to 3) 

9G jk 

3 y i =* G. ■ (i»l to 3 , j*l to 3) 

TE~ 13 

3 

4. 2. 7. 2 Attitude Transition and Covariance Propagation 

The gyro measurements can be used to calculate the state 
transitions and covariance propagations in a Kalman fil- 
ter. The attitude transition is given by 


q(t N+l>* ( cos( f At > m + - sin(^ At) [SI] } q(t N ) (4-38a) 
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= [Miq(V 
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(4-38b) 


Assuming that all the biases considered are constant in 
time and coupling between attitude and biases gives neg- 
ligible effects on state transition and covariance prop- 
agation, then the state transition matrix , the state 
noise matrix Q, and the state noise transition matrix 
? are identity matrices except for the three attitude 
components. Let q^ q 2 , q 3 be the three independent 
attitude components , then the corresponding <j> , Q , and 
v are given below. 
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(4-39) 

(4-40) 


where o* n • the gyro noise 

* a e + (T v At+ t c u (At) 3 

where a . 0 , and 0 are the standard devia- 

6 v n 

tions of the electronic, float torque, and 
torque derivative noise terius of the gyro. 
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(4-41) 
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4.3 OBSERVATION ERROR AND DATA MODELING 

The observation error matrix R in Eq. (3-1) and (3-2) 
should include the errors in the data due to both random 
noise and unmodeled systematic errors. However, the 
data errors due to unmodeled systematic errors are diffi- 
cult to estimate and require a great deal of computation 
time to calculate from point to point. Thus, in most 
practices, only the estimated random error is included 
in the matrix R. In this case, R is defined as a diag- 
onal matrix, with the diagonal elements equal to the 
square of the standard deviations of the estimated meas- 
surement errors. 

In the case of observability prediction studies, the 
actual or simulated data is not required by GCS. The 
only data information required by the predictor is a 
scheme to specify the times at which the data is avail- 
able. This scheme is referred to as "data modeling" 
method. For most practices, a simple data modeling 
method which assumes continuous data segments inter- 
leaved with specified data gaps is usually adequate to 
perform attitude and sensor bias observability studies. 

4.4 ATTITUDE VERIFICATION 

The attitude verification function compares the attitudes 
determined momentarily onboard with that determined from 
the ground using long data spans and recalibrated sensors. 
The standard deviation between the two types of attitudes 
is computed and compared with the ground attitude deter- 
mination accuracy as well as the attitude accuracy re- 
quirement. If this standard deviation is higher than 
both the ground attitude determination accuracy and the 
attitude accuracy requirement, then a flag will be set 


to indicate the need of updating the sensor calibration 
files to improve the onboard attitude determination ac- 
curacy . 


The standard deviation between the onboard deterministic 
attitudes and the ground attitude is obtained by the 
following equation M 


On = 


1 ( 9 i~ 6 g) 2 


N-l 




7a 


(4 


1 . 1 - 

where 0^ is the i pitch, roll, or yaw angle computed 
onboard, 6 g is the corresponding ground attitude, and 
N is the total number of onboard attitudes included in 
the computation. 


4-25 


- 42 ) 





REFERENCES 


1. Computer Sciences Corporation/ CSC/SD-78/6082 , Solar 
Maximum Mission (S'MM) Attitude System Functional 
Specifications and Requirements / R.W. Byrne/ P.J, 
Gambardella, J.A. Garber/ G. Nair, D.R. Sood, and 
G.K. Tandon, September 1978 

2. Computer Sciences Corporation, CSC/SD-79/6086, Solar 
Maximum Mission Attitude Determination System (SMM?ADS ) 
User*s Guide and System Description , J. Grondalski, 

V. Brown, H, Erickson, S. Hoven, R. Kwon, W. Lindboe, 

G. Nacios, and G. Neal, October 1979 

3. Computer Science Corporation, CSC/SD-79/6089 , Solar 
Maximum Mission Fine Pointing Sun Sensor (FPSS) off - 
Null Calibration System (SMM/FOCS) User's Guide and 
System Description , J~ Grondalski , J. Crowley, W. 
Lindboe, S. Watson, and D. Wilson, October 1979 

4. Computer Sciences Corporation, CSC/TM-80/6159 , Solar 
Maximum Mission (SMM) Attitude Analysis Postlaunch 
Report , G. Nair, R.H. Thompson, P.J. Gambardella, G.F. 
Neal, Y.R. Kwon, P.C. Kammeyer, and J.M. Buckley, 
August 1980 

5. James R. Wertz (ed.) Spacecraft Attitude Determination 

and Control, Boston: D. Reidel Publishing Company, 

1978 

6. Computer Sciences Corporation, CSC/TR-7 5/6001, Multi - 
satellite Attitude Determination/Optical Aspect Bias 
Determination (MSAD/DABIAS ) System Description and 
Operation Guide , Vol. 1, M. Joseph et.al., April 1975 


A-l 





